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\£) . We derive a set of exact cosmological solutions to the D = 4, M = 1 supergravity description of 

f*"^ ■ heterotic M-theory. Having identified a new and exact SU (3) Toda model solution, we then apply 

symmetry transformations to both this solution and to a previously known SU(2) Toda model, in 
CNl , order to derive two further sets of new cosmological solutions. In the symmetry-transformed SU (3) 

» I ■ Toda case we find an unusual "bouncing" motion for the M5 brane, such that this brane can be 

O |' made to reverse direction part way through its evolution. This bounce occurs purely through the 

, interaction of non-standard kinetic terms, as there are no explicit potentials in the action. We 

also present a perturbation calculation which demonstrates that, in a simple static limit, heterotic 
M-theory possesses a scale-invariant isocurvature mode. This mode persists in certain asymptotic 
limits of all the solutions we have derived, including the bouncing solution. 

CN ■ PACS numbers: 11.25.Mj, 11.25.Yb, 98.80.Cq 
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In the past ten years, heterotic M-theory has provided an exciting arena in which to analyse the cosmology and particle 
physics of our universe 

B H H 0, II Representing the low -energy limit of the strongly-coupled heterotic Eg x _E 8 
Oh, string, this theory not only combines gravitational, particle and braneworld physics into one unified description, but 
also possesses a detailed and constrained field-content that cannot be arbitrarily adjusted. Therefore, it retains a 
definitive and unambiguous relationship to M-theory itself. However, in the D = A,M = 1 supergravity description of 
heterotic M-theory, a number of important questions still remain unanswered. Consider, for example, the simple 
cosmological situation that occurs when we retain only the dilaton S, the universal T modulus, and the field 
Z describing a single M5 brane. Despite the fact that this leads to vanishing supcrpotcntial in four dimensions, 
the resulting cosmology is highly non-linear and demonstrates quite unexpected behaviour. In Ref. [3 the S, T, Z 
cosmology was analysed in a truncated limit with all axion fields removed. It was then shown that the scalar 
corresponding to the M5 brane position must be included in the set of cosmologically significant fields, and this leads 
to a forcing effect whereby the ambient dimensions change size as the brane moves. Moreover, the frictional forces 
acting back on the brane are such that it accelerates and then decelerates back to rest, mimicking a time-dependent 
force of finite duration. This illustrates the unconventional effect of non-standard kinetic terms in the theory, and 
the means by which the brane can undergo a single displacement by exchanging energy with its environment. In the 
special situation presented in Ref. this effect can be described exactly using the Toda formalism Q , and the 
model of Ref. is an SU(2) Toda model. 

Given this behaviour, it is interesting to consider whether more complicated trajectories for the M5 brane are possible. 
For example, not all of the scalar fields were considered in the SU(2) Toda model of Ref. Q, since the axionic fields 
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were consistently truncated away. This means that only a portion of the full solution-space was explored, and that 
the SU (2) behaviour is liable to be only an approximation once the axions are restored. Therefore, following on from 
the work of Ref. Q, we wish to analyse situations in which additional axionic fields are evolving in conjunction with 
the brane, and determine whether interesting new behaviours for the M5 brane can occur. In particular, we wish to 
determine whether an M5 brane can undergo multiple displacements, and even reverse direction in the absence of 
explicit potentials. 

Before embarking on the detailed calculations, we first summarise our results. We uncover a new and exact SU (3) 
Toda model, in which the M5 brane can undergo two successive displacements in the same direction. That is, the 
brane spontaneously accelerates twice in response to the other moduli fields to which it is coupled. Applying the 
symmetries derived in our companion paper Ref. to this model, as well as to the known SU(2) model of Ref. 0, 
we obtain two additional sets of new solutions. In the symmetry-transformed SU(3) case the brane can undergo 
two successive displacements in opposite directions, and so reverse direction and "bounce" without the presence of 
any explicit potentials in the action. This effect occurs purely through the interaction of non-standard kinetic terms, 
via the cross-couplings of the various fields, and constitutes an exact supergravity solution that has been rigorously 
deduced from M-theory. Finally we investigate the generation of density perturbations in these models, and show 
that heterotic M-theory possesses a scale- invariant isocurvature mode in some of the axion fields. This last result is 
consistent with the original findings of the pre Big Bang scenario 0, 0| and in agreement with the result obtained 
in Ref. |l3l | where it was first shown that the moving brane itself could not generate a scale-invariant perturbation 
spectrum. Following a conclusion, in an appendix we present the technical details of the SU(3) Toda model derivation. 



We now review the D = A,N = 1 supergravity action presented in Ref. pj- Recall that this was derived via a 
compactification of 11D supergravity on the orbifold S 1 /Z2 x CY3, where CI3 denotes a Calabi-Yau three- fold. This 
leads to two four-dimensional boundary planes separated along a fifth dimension. If the fifth dimension is labelled 
by a normalised coordinate z G [0, 1], then the boundaries reside at z = 0, 1 respectively and have the charges qo,qi- 
A single M5 brane is also included in the space, by wrapping it on a holomorphic 2-cycle of the CY3. The brane 
then appears as a three-brane of charge q that lies parallel to the boundaries, and which can move along the interval. 
Importantly, the interaction between the boundaries and brane leads to the existence of a static, triple-domain wall 
BPS solution. One can then consider further reducing on this solution, so as to find a supergravity theory describing 
slowly varying fluctuations about the static BPS vacuum. This contains the six scalar fields /3, x, 4>, °~i z i v with the 
following non-standard kinetic terms 



Each of these scalars has an underlying significance in terms of the D = 5 parent theory from which it descends. 
The scalar (3 is the zero-mode of the 555 component in the D — 5 metric, and measures the separation between 
the boundaries. Specifically, the separation is given by irpe 13 in terms of some dimensionful reference size up. The 
field <fi represents the orbifold-averaged Calabi-Yau volume, such that the physical size is given by ve^ in terms of a 
dimensionful reference volume v. The scalars er, x originate from the bulk three-form and graviphoton field respectively. 
The field z measures the position of the bulk brane between the boundaries, with the points z = 0, 1 corresponding 
to the boundaries. Lastly, the field v arises from the self-dual two-form on the brane worldvolume. 

This reduction on a BPS solution guarantees that the scalars must group into supersymmetric multiplets described 
by a supersymmetric action. One can verify that they naturally fall into the pairs (</>, tr),(/3, x)A z i v )i which are the 
bosonic components of chiral superfields S, T, Z as follows 



II. THE FOUR-DIMENSIONAL ACTION 




(1) 





Z = e fJ z + 2i(-v + zx) 



(2) 



This naturally leads to a Kahler manifold expression for the scalar part of the action 




(3) 
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where the superfields are grouped into a coordinate vector <I> = (5, T, Z), with the complex conjugate coordinates 
denoted by <E>. The Kahler metric Kq is given by 

d 2 K 



Ki3 



in terms of the Kahler potential 



K 



In 



, (Z + Z)^ 
b + b — q 



3 In (T + T) 



T + T 

This Kahler potential is computed only to linear-order in the two parameters (fc = 1, 2) defined by 



(4) 



(5) 
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Here k is the eleven-dimensional Newton constant, and irp, v are the dimensionful scales mentioned above. The two 
conditions eu <C 1 then restrict the accessible regions of moduli-space in which we can trust the four-dimensional 
effective theory. In addition, the supcrgravity action Eq. ^ can only be trusted in the limit where stringy a' 
corrections arc suitably small, as these corrections introduce higher-derivative terms that we have disregarded. 



III. EQUATIONS OF MOTION 

We now turn to the equations of the motion arising from the action Eq. Q . If we assume a spatially flat Friedmann 
Robertson Walker (FRW) cosmology for the four-dimensional spacetime, then the metric takes the form 

ds\ = ~e 2n{T UT 2 + e 2a ^ T) d ij dx i dx j (6) 

where i,j = 1, ... ,3, the scale-factor is a(r), and n(r) represents a gauge freedom in the choice of time coordinate. 
Denoting a r derivative by an overdot, and assuming all fields are purely functions of r, one obtains the Einstein field 
equations 

-3a 2 + + - A (i 2 + \qe^z 2 + Se^ X 2 + \e~ 2 * (a + Aqzvf + 2qe~^(y - X z) 2 = (7) 
2d + (3d - 2h)a + + -Ji 2 + ^qe^z 2 + 3e~ 2t3 X 2 + \e- 2 ' i ' (& + AqzO) 2 + 2ge"' 3 ~^(i> - X zf = (8) 

the 0, /9, z equations of motion 

4> + (3d - n)4> + qeP'^z 2 + (a + Aqzvfer 24, + Aq{v - xzfe' -' 1 ' = (9) 
3/3 + 3(3d - n)/3 - qe^z 2 + 12xV 2/3 + Aq{v - xzfe~ ~' l> = (10) 

J- { [e^z - 4 X (z> - X z)e- f3 -' t >] e 3a - n } e™" 3 " - 2v{a + Aqz^e- 2 ^ = (11) 
and the a, v, X equations of motion 

4- { [ (& + ±qzv)e- 2 * 1 e 3a ~ n } = (12) 
dr 

4 { \z(a + 4qzv)e- 2 * + 2{v - X^e'^] e 3 "-"} = (13) 
dr 

4 [ (3e~ 2l3 x)e 3a ~ n 1 e™" 3 " + 2qz(v - X i)e~' 9 -* = (14) 
dr 

As we cannot exactly solve these equations of motion, we will utilise the following two solution methods. Firstly, we 
search for specialised solutions that occur when the equations are truncated, usually by setting certain combinations of 
fields to zero. This will allow us to recover the known SU(2) Toda model of Ref. Q, as well as a previously undiscovered 
SU(3) Toda model. Secondly, we will utilise the scalar-field symmetry transformations that were derived in our recent 
companion paper Ref. |10( . That is, we will apply these symmetry transformations to the fields of the SU(2) and 
SU(3) Toda models in turn, and so derive two new solutions to the equations of motion. We will find that in these 
new solutions the M5 brane can evolve in far more complicated ways than has previously been seen. 
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IV. REVIEW OF THE SU(2) TODA MODEL 



We now briefly recall the behaviour of the SU{2) Toda model found in Ref. 0. This will prove useful because the 
SU(2) model exhibits features that persist in all the solutions we will present, and so will illuminate the discussions 
to come. In addition, it is worthwhile studying this model in order to understand how symmetry-transformations will 
affect it. 



A. The S*7( 2 ) Toda Model Solutions 



The SU (2) model can be derived by choosing the axions to satisfy a, v = constant, % = 0. The remaining fields 
a, <f>, j3, z can then be solved for exactly, essentially due to the fact that the field z satisfies the conservation law 



/3-ip+3a-n ■ _ 



constant 



(15) 



Namely, inserting this result back into the remaining equations of motion yields a closed set of equations in a, j3, 4> 
which can be solved in isolation. In particular, these equations can be reformulated in terms of the motion of a 
"particle" with coordinates a = (a, (3, <ft) which roams over a three-dimensional space and experiences an exponential 
potential U . These equations take the form 



4- {EGa) 

CLT 



oa 



1 T 

-EotGa 

2 



E~ X U = 



which can be derived by variation of a, E from the simpler particle Lagrangian 

£ = -Ea T Ga - E~ X U 



(16) 



(17) 



Here we have defined a moduli-space metric G = diag(— 3, |, j) of Minkowski signature, and a particlc-worldlinc 
metric E = e -' l+d a where d = (3, 0, 0) is a dimension vector that gives the number of spatial dimensions associated 
to the scale-factor a. The function E thus encodes the arbitrary choice of time parameterisation of the worldlinc, 
with its variation in the Lagrangian naturally producing an energy conservation constraint. Finally, the moduli-space 
potential U = U± is given by 



I exp (qi 



a ) 



qi = (0,-1,1) 



(18) 



where u\ is a positive constant. Up to a constant length rcscaling, the vector qi defines the single, simple root-vector 
of the Lie algebra SU(2), and so the Lagrangian Eq. i|17|) defines an SU(2) Toda model. This is an exactly intcgrablc 
system, and in the proper-time gauge n = one finds the general solution: 



a — oto = pi In 



t-t 
T 



(Pf - Pi) In 1 





t-t 


1 




T 





(19) 



zo 



d 1 



t-tp 
T 



(20) 



p 7 Gp 7 = 



p 7 d 



qi • a = In 



qd 2 < qi , qi > 



S = qi • Pi 



Pf ~Pi = 5 



< qi,qi > 



(21) 



Here the subscript 7 takes the values 7 = i, f, and the scalar product < •, • > is defined by < a, b >= a T G _1 b. Note 
that the constraints in Eq. (|21|) must be enforced so that Eqs. i|19fl - H2(J|) are indeed the correct field solutions. Once 
this is done, the solution describes a transition between two asymptotically free-field states. That is, the initial field 
velocities are equal to the "expansion power" constants Pi = (1/3,^^,^^^), the final field velocities are equal to the 
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constants Pf — (1/3, Pp,f,p<j>,f), and the nonsupersymmetric second term in Eq. I|19|) forces a smooth acceleration 
between these "rolling-radii" (rr) regimes. In fact, the underlying reason for this Pi — » p/ interpolation is the motion 
of the brane itself, which according to Eq. H2()(l is at rest in the extreme limits, but borrows kinetic energy from /?, <j> 
and moves significantly at the intermediate time t — to « T. For the sake of concreteness, we now present the explicit 
form of the solutions by inserting the various vector quantities. This gives 



a ~ a ° = 3 ln 



t-t 



f3 - /3o = Pp,i In 



In 



T 

t-to 



T 

t-to 



T 





t-to 




1 




T 




(» 


t-to 


1 




T 





-1/5 



These are subject to the relations 



: Pf3,i - P<f>,i 



as well as the "ellipse" constraint 



In 



2q<P 



P0J 
P<t>J 



4 

3 



1 1 

3 -1 



P/3,i 
P<t>,i 



(22) 



Notice, in fact, that if we enforce the constraint that fixes (Pf3,f-,Pc/>,f) in terms of (pp.i,P<t>.i) then the final expansion 
powers (p/3j,P<f,j) are automatically guaranteed to lie on the same ellipse. The interpretation of this ellipse condition 
is relatively simple, and can be illustrated in a phase-plane diagram as follows. Consider drawing a 2D plot where the 
horizontal axis corresponds to d[3/du (with u = 3a) and the vertical axis to d<p/du. Then the constants (p/3, 7 , J30, 7 ), 
where y = i, /, are the asymptotic values of d(3 / du and d(f>/du. Thus, they correspond to the values attained in the 
plane at the extreme endpoints of the phase-plane trajectory. Hence, according to Eq. (|22|l and the constraints, the 
trajectory traced out in the phase-plane must begin and end at two different points on a single, fixed ellipse drawn in 
that plane. These phase-plane plots, or "ellipse diagrams" as we shall call them, will prove to be extremely useful in 
exhibiting the behaviour of the system diagrammatically. This is because the shape and curves of the trajectories in 
these diagrams tell us very graphically about the brane motion and changes to the axions. 



B. Analysis and Validity of the 5(7(2) Toda Model 



To exhibit the behaviour of the fields, we now plot an ellipse diagram. This proves to be far more intuitive and 
useful than following the behaviour of all fields individually. Before doing this, we must recognise that the solutions 
in Eq. (|19|l are valid over two disconnected time ranges given by 

i c / t -00 '^) ' (-) branch , . 

te \ (i ,+oo) , (+) branch [Z6) 

where the time t = to corresponds to a curvature singularity. Consequently, there are two different notions of "early" 
and "late" built into the solutions, depending on the choice of branch. For example, although t — to corresponds to 
a past singularity in the (+) branch, it corresponds to a future singularity from the perspective of the (— ) branch. 
The (— ) branch is, in fact, an example of a pre Big Bang (PBB) era, which is automatically undergoing supcrluminal 
deflation. 

Therefore, to avoid confusion we must always pick a particular branch, and take care with what constitutes early and 
late behaviour. In particular, the Pi constants only correspond to an "initial" set of expansion powers as implied by 
the subscript if they satisfy 

r f > (— ) branch ,„ 

d= ^~^M<0 (+) branch ^ 
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That is, only those powers satisfying this condition can ever be early time states of the system. In Fig. E]we have 
plotted some representative trajectories on the (— ) branch. Note that the fields /3,<f> start at a single point on the 
ellipse, with their initial powers corresponding to that sector with S > 0. On the negative branch this "early" time 
state corresponds to the infinitely negative past t — to — > — oo. The fields then evolve such that the effective trajectory 
in the plane is a straight line. This linear behaviour is a consequence of the relation 3df3/du + dcj>/du oc constant, and 
this in turn is possible because all the axions have been truncated away. The trajectory then ends on the opposite 
sector of the ellipse, ending up in a "late" time state as t — to — > from below. This directed evolution between parts 
of the ellipse cannot be reversed unless we switch the branch from (— ) to (+), so the accessible early-time states of the 
system are fixed by the choice of the branch. Thus, interesting physical results will sometimes necessitate choosing 
one branch over another. 




8>0 



Pp 

FIG. 1: Figure to show the linear mapping across the ellipse, with the direction fixed by the choice of time branch. If we pick a 
candidate set of expansion powers satisfying S > 0, then these are indeed available early-time states on the (— ) branch. States 
satisfying 8 < are not available at early time; instead, they are the late time states that the system evolves into. This is all 
reversed on the (+) branch. 



To complete this section, we now comment on the validity of these SU(2) Toda solutions. Recall that the four- 
dimensional action Eq. Q is known only as a power series, with five-dimensional gravitational corrections measured 
in powers of the (k = 1, 2). In the above model one finds that 

J \t — to\ +s — ► oo early time 
tk ~ I \t~to\- & -> oo late time [ ' 

This divergence is due to the fact that the coupling of the bulk brane to f3, <j> is itself proportional to the e^, so that 
when the brane moves the system is necessarily driven to a five-dimensional limit in both asymptotic regimes. This 
will cause the four-dimensional theory to break down, and with it the solutions Eq. (|19() . This divergence is in fact 
familiar from PBB cosmology where the addition of the dilatonic axion to the dilaton-moduli system causes the same 
problem. A more insidious problem, however, is that even at intermediate times one cannot make <C 1 whilst 
simultaneously fitting the entire z displacement profile within the physical orbifold extent z £ [0, 1]. Say, for example, 
that we search for the minimum value of the parameters. One can verify that this occurs at precisely the time 
t — to = T, and at this time the ellipse trajectory intersects the line 6 = 0. The magnitude at the minimum is then 

^fc|min ^ ,9 



Evidently, d 3> 1 is now required in order to make this small. Having made this choice, there will be a finite period 
of time, depending on the value of d, where the are small and our solutions, Eq. (|19|) . arc valid. This period 
ends with the rapid collision of the brane with the boundary which, of course, also invalidates the above analytical 
solutions. (For a discussion of the evolution after the collision see Ref. [13 )■ The limitations on the validity of these 
solutions, both due to the constraint and brane-boundary collision, may be viewed as a disadvantage and one may 
ask whether other solutions with a larger range of validity exist. We will show that this is indeed the case for some 
of the new solutions to be discussed below. 
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V. THE SU(3) TODA MODEL 



After reviewing the SU (2) model at some length, we will now present an entirely new SU (3) solution. This proves 
to be significantly more complicated than the SU(2) Toda model, as wc might expect from the fact that the Lie 
group SU(3) is more complicated in structure than SU(2). We will find that the 5*7(3) solutions are fundamentally 
controlled by two influences, one due to the motion of the M5 brane and the other due to changes in \. This leads 
to new, characteristic trajectories in the ellipse diagrams. In particular, the coupling between z and \ allows the M5 
brane to undergo two successive displacements in the same direction. This generalises the behaviour of the old SU (2) 
case, and demonstrates that the brane can undergo repetitious displacement behaviour. 



A. The 5*7(3) Toda Model Solutions 



To derive the SU (3) Toda model, one chooses the axions to satisfy 

v - xz = a + Aqzv = (26) 

This effectively sets two of the terms in the action to zero for all time, without forcing any of the fields or their 
derivatives to be individually zero. To see that this corresponds to an ,517(3) Toda model we impose Eq. (|26[) . and 
note that the z and \ equations are now total derivatives 

4- [e^ +3Q -"il = 4- \e- 2 P +3a ~ n x] = 
or dr 1 J 

These can be immediately integrated to give constants of the motion. Inserting these conservation laws back into 
the remaining equations of motion then yields a closed set of equations in a, (3, </> that can again be derived from the 
particle Lagrangian 

C=^Ea T Ga- E~ X U (27) 

The quantities a, G, E remain unchanged from the SU(2) Toda model, but the potential is modified to U = U\ + U2, 
with 

U\ = \^u\ cxp (qi ■ a) , qi = (0,-1,1) 

U2 = exp (q 2 ■ a) , q 2 = (0,2,0) 

This means that the effective particle motion of a is now subjected to two exponential forces. To be a Toda model, 
a precise relationship must exist between the orientations and lengths of the vectors defined by the q^. Consider the 
following matrix 

< q*,qj >= 3 ' ( -1 2 ) (28) 

where the scalar product is once again defined by < a, b >= a T G _1 b. The right hand side of Eq. I|28|l happens to 
be a constant multiple of the Cartan matrix of the SU (3) Lie algebra, so that up to a constant length rescaling the 
vectors qi, q 2 are identical to the two simple root-vectors of 5**7(3). Thus, the model is an exactly integrable SU(3) 
Toda model. In particular, an exact, analytical description of the behaviour is now accessible if we choose a basis for 
the moduli-space that is adapted to the SU (3) root- vectors. This decouples the equations of motion and allows them 
to be readily solved, the complicated details of which are reserved for Appendix A. The proper-time solutions for the 
fields in the gauge n = are then given by 



a — otQ = pi In 



t -tp 
T 



(P/ ~Pi) ln 



t - t ( 
T 



-6 



t-k 



-s 3 



t-t 




t-to 


T 




T 



-<5N 



-l/<5<9 



(29) 
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t-to 
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el 


t-to 


-Vv _1 


(» 








T 

















X - Xo = d x 1 + 



t-to 



T, 







t-to 



T 



-5> 



• 1 



t-to 



T 



(30) 
(31) 



The constants are subject to the following two sets of "5C/(2)-like" constraints 

'qd 2 < qi,qi > 



p 7 Gp 7 = p 7 • d = 1 S = — qi Pi qi • ccq = In 



8 



pf ] Gp { f x) =0 pf -d = l <5/3 = -q 2 -p t q 2 



a - Pj In 



( 34 < q 2 ,q2 > 



where 7 = i, /, and the scalar product < •, • > is again defined by < a, b >= a T G 1 h. Moreover, p/, p^and p^ are 
related by the two SU(2) maps 



(32) 
(33) 



2G-iqi (x) 2G" 1 q 2 



< qi,qi > 



< q 2 ,q 2 > 



(34) 



Finally, the fractional quantities 8 Z , 9 X are fixed according to 

qi ■ Pi 



q 2 ■ Pt 



(qi + q 2 )-Pi x (qi+q 2 )-Pi 



(35) 



These satisfy < 9 Z , 9 X < 1 and Z + 6 X = 1 . For clarity, we now present the solutions and constraints for a in 
component fields. These read 



a — oto = — In 



(3 - 00= P0,i In 



t-to 



T 
t-to 



T 



+ (PPJ ~ P0,i) ln 



+ (pj$ " P/9,i) In 



- 4>o = V<i>,i In 



t-tr 



T 



+ (P^'f - i>*,0 In 



t-to 


(i + el 


t-h 


T 







t-h 



T 







t-tr 



t-h 



T 



T 

t-to 



-6 



i + oi 



t - 1 







Sg 



T0 
t-to 



T 



-1/6 



-\-ys 



-1/8 



(36) 



(37) 
(38) 
(39) 



The constants 6,ppj,tQ,T,d, z all occurred in the previous SU(2) solutions and so are familiar. The three new 
constants are given by Tjg, xo> d x with the remainder constrained according to 



Sp = -2pp. 



I ) w Z 



5 + 5 f3 



Po = ln(2d x ) +pp,iki 



T 
T~0 



A) - 00 = In 



2qd 2 



P0J 
P<t>J 



1/1 1 
2 V 3 -1 



P0,i 



1 
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Before discussing these SU(3) solutions in more detail, we should also comment on the solutions for the additional 
fields v, a satisfying Eq. I|2()|) . It transpires that the a solution involves a non-elementary integral, and so cannot be 
presented analytically. However, its behaviour can always be computed numerically. On the other hand, the field v 
takes the simple form 



t-t 



T 



t-t 



-5 b 



t-t 



-5 B 



(40) 



subject to the conditions 



Xo (xo 



(41) 



Interestingly, the field v can reverse field- velocity midway through its evolution, and so "turn-around" or "bounce" . 
This peculiar effect will have important ramifications when we transform the SU (3) solutions later on, for it will allow 
the branc field z to bounce as well. 



B. Analysis and Validity of the SU(S) Toda Model 



There are two distinct SU(2) models embedded non-trivially in these solutions. If we take the limit Tp ^> T then the 
early-time behaviour of the fields is formally identical to the solutions Eq. (|19fl of the previous section. If instead we 
reverse the temporal sequence by choosing T^«T then the early time behaviour is another three-field SU (2) model 
involving \ as the axion. These two models are not decoupled as they would be in the SU(2) x SU{2) Toda case, but 
instead are non-trivially mixed inside the SU(3) model. Only in extreme cases can we discern the underlying SU(2) 
components, and so at a general, intermediate time there will not be a clean separation of the effects of the z and 
X motions. Indeed, these two embedded behaviours couple and compete with one another, and attempt to drive the 
expansion powers according to the two conflicting processes 

M 

Pi -> Pf , Pt^ P f 

In general, this means that neither p f nor pi individually succeed in becoming the actual expansion powers that the 
system adopts at late-time. However, the system may spend some part of its evolution at intermediate rolling-radii 
states where it adopts these powers temporarily. The true late-time rolling-radii states p'^ are in fact determined from 
a "combined" relation given by 

2G _1 q 

p' f - Pi = (6 + Sp) , q = qi+q 2 (42) 

< q, q > 

Notice that the final states are computed as if the system followed an ordinary SU(2) model with a combined 
parameter 5 + Sp. However, the intermediate behaviour strongly deviates from any such simple SU(2) evolution, and 
we should treat Eq. (|42|l merely as a formal tool for deducing the rolling-radii endpoints of the trajectory. 

To see this clearly, we plot the field behaviour on the ellipse as we did in Section llVl (see Fig. 0). Again, the solutions 
break in (±) branches, with p^ corresponding to the early-time expansion powers only if they simultaneously satisfy 
the two inequalities 

8,5$ >0 on (-) , 5,5 P <Q on (+) (43) 

The additional 8p condition restricts the accessible early-time powers to a narrower range of states compared to the 
SU(2) model. In Fig. [3] we also plot the displacements of z and the x axion. This illustrates the important fact that 
the fields z,x can undergo two successive displacements, since each is coupled to the time-development of the other. 
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FIG. 2: A typical set of SU(3) trajectories. Curve A has two horizontal lines and one diagonal line, and so contains two \ 
displacements and one z displacement. Curve B represents a special, degenerate case for which both z and \ evolve at once and 
mimic a single field. Curve C has two diagonal lines and one horizontal line, and so contains two z motions separated by a x 
displacement. All intermediate cases between these curves are possible. 




FIG. 3: An interesting example of the field evolution. The left hand plot shows the fields (5 and <j>. The kinks in these curves 
are caused by the displacements of z and x, which are shown in the right hand plot. Notice that z moves twice in succession. 



To complete this section, we now comment on the validity of these SU (3) Toda solutions. One can show that 
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t-tp 

T 



(44) 



Using this, one can easily verify that in the asymptotic limits 



t-tp 
T 



+6 



oo at early time 



efe 



t-to 



oo at late time 



Notice that this follows only because 8, 8p are both positive (negative) on the negative (positive) time-branch. 
Consequently, the SU(3) solutions cannot be trusted asymptotically, as with the previous SU(2) solutions of 
Section llVl Further investigation of Eq. (|44|l also reveals that it can never be made smaller than the leading coefficient, 
which is of order 1/qcP. This demonstrates that the smallest attainable value of the ek is given by 



efc|min~-4- (45) 
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Hence, to achieve <C 1 the solutions require us to take d^$> 1 and allow the brane to leave the orbifold interval. As 
such, the SU (3) model has similar problems to the SU (2) model. Of course, as long as we are interested in relatively 
short timcscales, and always concentrate on the brane behaviour inside of the interval but away from the boundaries 
(and the collision) then no particular problem is posed. Away from the boundaries the SU(3) solutions with d 3> 1 
are reliable for a short time, and the fact that the brane must eventually leave the interval does not change this fact. 
Hence, there are always regions where all fields are evolving in an eu *C 1 regime with the brane inside the interval. 

However, it would obviously be valuable if these regions could be extended to cover the entire displacement profile of 
the brane, such that the brane moves and comes to rest while remaining inside the interval with < 1 throughout. 
Although this is impossible with the SU (3) solutions themselves, when we come to symmetry-transform the SU (3) 
solutions we will find circumstances under which d < 1 and e& <C 1 simultaneously. 



VI. APPLICATION OF THE SYMMETRIES 



We now apply the symmetries presented in our companion paper Ref. jlCj to the 577(2) and SU(3) models in turn. (For 
previous work in this area, also see Ref. |l5l|\ These symmetries mix the scalar fields together in new combinations, 
and yet leave the action Eq. Q invariant. Consequently, the new time-dependent combinations for the fields that 
emerge, no matter how complicated, still solve the equations of motion. The seven-dimensional symmetry group G is 
a maximal parabolic subgroup of Sp(4, R) given by 

G = SL(2, R) k SU(1, 2)/U{2) (46) 

where k denotes a scmidircct product. The scalar-field space is then the group manifold M. = 5p(4, R)/C/(2), but 
equipped with a nonhomogencous Ricmannian metric such that the total symmetry group fills out only G and not 
the whole of 5p(4, R). There are seven, distinct types of transformations L l that can be applied to the scalar fields, 
each one controlled by a continuously adjustable real constant c^. These are given by 



L 1 : 
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X^ — 7Z-, so , i 9 9,3 o^o-c 2 -2qv z^z + c 2 v 



\(*eW (l + c 2X ) 2 + icle^ 

a — > a — C3 • 2qz 2 v — > ^ + C3Z 



z -> z e C4/2 1/ -> 1/ e C4/2 



Z + C 5 



(47) 

Some of these transformations will simply amount to rcparamctcrisations of the existing integration constants in the 
solutions they are applied to. Some, however, change the solutions into new functional forms, which will in turn be 
completely new solutions to the equations of motion. In particular, note that L 1 , L 2 , L 3 represent the SL(2,M) group 
of transformations, which act on S, T, Z as follows 

T ' =C §F^1 ' Z ' = idFTd ' S ' = S -J§Tl With a > 6 ' c ' deR > ad-bc=l (48) 

Here the four constants a,b,c,d (subject to one constraint) are proportional to c\, 02,03, and are better adapted to 
the SL(2,M) symmetry In particular, they can be considered the four entries of a 2 x 2 SL(2, R) matrix. Having now 
rewritten the SL(2, R) transformations in this compact form, note the crucial fact that Eq. I|48() does not represent 
ordinary T-duality, for the complex coordinates S and Z must also be transformed. Nonetheless, the action on T alone 
is indistinguishable from conventional T-duality, and so all-told we will dub this a "generalised" T-duality. These 
generalised T-duality transformations can produce exceedingly complicated new behaviours, and can significantly 
affect any existing time-dependent solutions that they are applied to. Consequently, we can expect to derive new 
solutions to the equations of motion by transforming the SU(2) and SU(3) models using these symmetries. While 
it was shown in Ref. |lfjj that the symmetries do not form a transitive group on Ai, so that we cannot use them to 
build the general solution to the equations of motion, we can nonetheless make significant progress in this direction. 
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VII. TRANSFORMING THE SU{2) MODEL 

We now apply the finite symmetries L % to the SU(2) model. Only L 2 ,L 3 have a non-trivial effect, and they modify 
the system such that it is no longer a simple Toda model that can be solved using the Toda methodology. This, of 
course, is the crucial reason why we use the symmetries in the first place, as they allow us access to complicated 
new solutions that we cannot otherwise uncover using standard methods. Although the brane only undergoes one 
displacement, we will find that the parameters can have significantly different development in these transformed 
solutions. Specifically, in certain cases the ek are naturally decreasing into the past or future. 



A. Transformed-5'L r (2) Solutions 

One can verify that the symmetries transform the SU (2) solutions into the following form 
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(49) 
(50) 
(51) 

(52) 
(53) 



v — uq = 

Here we have defined the combinations 

s = ±l , 2A = l-s , 5 = 5 + A5p 

We are, however, not free to pick s in an arbitrary fashion, as the choice of sign crucially depends on the choice of 
initial expansion powers. The permissible choices are listed in tabled 



Expansion power 
range 


8,6p < on (+) 
S, 8,3 > on (-) 


5 + 5/3 < 0, dp < 0, 8 > on (+) 
5 + 5 > 0, dp > 0, 6 < on (-) 


Allowed choices 


s = ±1 


s = -1 



TABLE I: Table to show which of the values s = ±1 are valid choices given a set of initial expansion powers. Note that s = +1 
is not a valid choice in the second row of the table. 



The remaining constants are then subject to the same constraints as in the SU(3) model, namely 

'qd 2 < qi,qi > 
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W-d = l 



5/3 = -q2 • Pi 
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3dl < q 2 ,q2 > 



(54) 
(55) 
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where 7 = i, f and p/, p^and pt are related by the two SU(2) maps 



D (x) 



2G~ 1 q 2 
< q2,q2 > 



P/ 



^^q_ 
< q,q > 



q = qi + Aq 2 



(56) 



Notice the crucial fact that the system does not necessarily have the same asymptotic behaviour as the SU{2) model. 
To see this, we note that the early-time powers pi can now be taken from anywhere in the region 

5 + 6p,5p<0 on (+) , 5 + 6p,5p>0 on (-) 

Most of this region was unavailable in the original SU(2) model, and so the permissible asymptotic behaviours have 
expanded into completely new regions. This is very different to the SU (3) model, which consistently narrowed the 
range of powers compared to the old 5(7(2) case, but did not expand the allowed range of powers at all. These newly 
accessible regions are entirely a consequence of the symmetry transformations, whose effects were entirely absent in 
the original SU (2) and SU (3) cases. Therefore, we can anticipate entirely different behaviour for the parameters 
in the asymptotic limits. 



For clarity, we now present the component field representation of a: 
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(57) 
(58) 



(59) 





As in all the cases considered, the expansion powers (pa,i,P<j>,i) are constrained to the ellipse as defined in Eq. (|22|l . In 
combination with the constraints above, this automatically forces /3 and <p to be in rolling-radii regimes at late-time 
that are also on the ellipse. 



B. Analysis and validity of the transformed-S'C/(2) model 



As in the previous section, we plot a particular example of the /3, 4> evolution across the ellipse on the (— ) branch 
(see Fig0J. The behaviour breaks down into three generic cases based on the relative magnitudes of the timescales 
T and Tp. Notice that, irrespective of these magnitudes, the field z can only ever undergo one displacement, and so 
behaves in a manner identical to the old SU (2) case. However, the crucial thing is that we can now achieve the same 
SU(2) behaviour for z inside a set of solutions that have completely different development for the In the particular 
example given, the originally diverging values of ek at late-time are now decreasing to arbitrarily small values into 
the future. Thus, the solutions become more and more reliable into the future. This is in stark contrast to the SU(2) 
model from which they originated, and demonstrates that the new \ behaviour is crucial in suppressing gravitational 
corrections to the four-dimensional theory. 
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FIG. 4: Transformed 5(7(2) ellipse behaviour for Tp T , Tp ~ T, Tp <C T, plotted on the ( — ) branch. The trajectories begin 
at the lower left as t — to — > — oo, and evolve to the upper right as t — to — > 0. 



Moreover, in Fig. |SJwe plot the Tp <C T field behaviour of z,x, such that the z displacement occurs in the reliable 
e/c <C 1 regime after the change in \- Crucially, this means that the brane motion can occur in a < 1 region without 
requiring d > 1. Such behaviour could never have occurred in the original SU(2) model, and is a consequence of the 
manner in which the fields z and x are incorporated together into a new, global structure for the overall solutions. 




IT/tl IT/tl 

FIG. 5: In the left hand plot we see <j> an< l the SU (2) -transformed solution for (3, with the remaining curve corresponding to the 
original untransformed SU(2) solution for /3. In the right hand plot we see the distinct z and \ displacements, with the motion 
of the bulk brane z occurring second. This allows the brane to displace with d < 1 and yet still be in a 6k >C 1 regime. 



The fact that the displacement of z can be made to entirely occur in a C 1 regime, without requiring d > 1, 
constitutes a significant improvement over the original SU{2) model. This proves that a reliable solution for z does 
not necessarily require it to eventually leave the compact space . However, there are obviously other possible examples 
beyond those shown in Figs. l4"I51 and we should now clarify the precise circumstances in which the can be made to 
decrease. Once again we consider the functional form of exp(/3 — cf>): 
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(60) 



By choosing signs appropriately, either the early or the late time limit can become "weakly-coupled" with etCl. 
However, only one of the asymptotic limits can be weakly-coupled, with the other still becoming "strongly-coupled" 
with efc S> 1. There are, of course, still solutions where efe ^> 1 is attained in both limits. The full state of affairs is 
summarised in table ITTl 
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s = +1 


s = -1 


s = -1 


o, Op < on (+) 
5,(5/3 > on (-) 


0, Op < on (+) 
5,(5/3 > on (-) 


o + Op < 0, dp < 0, d > on (+) 
<5 + (5/3 > 0, <5,g > 0, 6 < on (-) 


SW 


SS 


ws 



TABLE II: Table to show the asymptotic values of the coupling parameters et, depending on the sign of s. The notation is as 
follows: strong-strong (SS), strong-weak (SW) and weak-strong (WS), where the first word corresponds to the values in the 
early-time limit, and the second refers to their values in the late-time limit. 



Thus, there are three distinct types of solutions: weak-strong (WS), strong-weak (SW) and strong-strong (SS). In all 
three cases we can arrange for the z motion to occur in a £/; < 1 region with d < 1. To do this, we simply recognise 
that at t — to w T we can always take T@ <C T, and this will decrease the values of the ek below 1 without requiring 
ti> 1. Consequently, there is a tremendous degree of flexibility in the solutions, and cases with ei « 1 and d < 1 
are quite generic. 

Before leaving this section, we should also comment on stringy a' corrections. These become strong as we probe small 
length scales at (3 ~ 0, and so encounter new physics not accounted for in the effective supcrgravity description. As 
such, one must always ensure that (3 ^S> to trust any supergravity solution. We note that this is always possible for 
certain periods of time by an appropriate choice of integration constants, and so there is no obstruction to finding 
regimes where tk "C 1 and a' corrections are extremely small. The transformed SU(2) solutions thus incorporate all 
of the z behaviour from the SU (2) model, but now allow it to be compressed inside of the orbifold interval whilst 
simultaneously suppressing all unwanted corrections. 



VIII. TRANSFORMING THE SU(3) TODA MODEL 

We now apply the symmetries L l to the SU(2>) solutions, and so find a further class of new solutions. One finds in 
this case that only the action of the L 2 transformation can ever lead to new behaviour. This can be easily understood 
by noting that all the other transformations leave the SU(3) truncation conditions Eq. (|26fl invariant, while L 2 allows 
v — \z to become non-zero. The "activation" of this combination takes us outside of the original SU (3) Toda model, 
and into a new situation that is not itself solvable by Toda methods. Nonetheless, the symmetries allow us to access 
an exact, analytical description of the behaviour when this combination is non-zero. We will find that the brane 
can undergo two displacements in opposite directions, and so reverse direction without the presence of any explicit 
potentials. We will often call this a "bouncing" solution. 



A. Transformed-S'l7(3) Solutions 

These new solutions, although exact, are complicated and difficult to present in an elegant fashion. One means of 
presentation is to utilise two time-dependent functions p, r that are implicitly defined via the relations 

M4 + P 2 )- 1 =e \su {3 ) , priA+p 2 )- 1 = X \su(3) (61) 

These are built out of the (3, x solutions from the old (untransformed) SU(3) model. The new transformed-5f/(3) 
solutions can then be written in the form 

a = a|s(7(3) , <f) = (f>\su(3) , V = v\su(3) 

-1 

, z = z\su(3) + c 2 v\su(3) (62) 



(3 = In ^ 4r 



4 + (p + c 2 r)' 



X = r(p + c 2 r) 



4+ {p + c 2 rY 
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Here Ci is the real constant associated to the L 2 symmetry, and so corresponds to a new integration constant that can 
be varied at will. The remaining constants, it should be emphasised, are taken from the original ST/(3) model, and 
we should treat their values as determining an embedding of the old SU (3) behaviour inside the newly transformed 
solutions. Indeed, the fields a, 0, v are unaffected by the transformations, and evolve as in the old SU(3) case in any 
event. 

Notice as well that we cannot present a analytically, due to the fact that the corresponding SU(3) solution can only 
be computed numerically. As such, the transformed a solution must also be computed numerically. However, we 
emphasise that these numerical computations can be readily carried out with no obstruction, and that the symmetry 
transformations induce perfectly sensible behaviour for a in all cases. In addition, a can have no bearing on the 
time-development of the other fields, as the condition a + Aqzi> = is preserved under the symmetry group Eq. l|46f> . 
This means that a never appears in the equations of motion of the other fields, and so can never induce any changes 
to the brane or remaining axions. Consequently, we will not particularly concern ourselves with a from this point on. 

Having derived these new solutions by applying the symmetries, we must now consider the ramifications for the 
various fields involved. In particular, we are most interested in the field z and the issue of whether we can now achieve 
a sensible displacement at weak-coupling. In the original SU (3) model the field z could undergo two successive 
displacements in the same direction, but it could not do so whilst entirely within a ft « 1 regime. In the above case, 
however, the new field z is an additive mixture of the old SU(3) behaviours for z,v. This creates a significant new 
level of flexibility, and in the next section we will investigate the consequences for the brane and its displacements. 



B. Analysis and Validity of the Transformed-S'(7(3) model 

Due to the complexity of the solutions, the field behaviour is somewhat difficult to determine by mere inspection. 
However, one can verify that the symmetry-transformation does not affect the asymptotic development, and so the 
same set of states are accessed on the ellipse at early and late time as in the SU(3) model. However, the intermediate 
evolution is substantially, and interestingly, different. In Fig. we plot some examples on the ellipse. 
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FIG. 6: Here we show some examples of the transformed- SU (3) solutions. Curve A is the original SU(3) solution, and this 
can be progressively shifted toward curve B and into curve C as we change the integration constants. 

The particularly interesting feature of these new solutions is the motion of the brane. Specifically, for certain special 
choices of constants, the brane can "bounce" and spontaneously reverse direction midway through its evolution. 
Moreover, a thorough investigation of the parameter space reveals that it is possible to make <C 1 whilst z is 
undergoing this bounce strictly inside the orbifold interval. This is shown in Fig. [7| 

To see that this behaviour is indeed a consequence of the solutions in Eq. ijoTJl. one can proceed in the following 
qualitative fashion. First, we recognise that the effect of the symmetry transformation is to switch-on the combination 
v — xz to a non-zero value. This then acts as a driving force that modifies the original SU(3) evolution of z. Secondly, 
we recognise that this "modification", at a practical level, amounts to additivcly mixing the SU(3) behaviours of z 
and v together (see Eq. (03). So not only is z affected, but it is affected by a non-trivial mixing together with the 
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FIG. 7: 27ie piot shows the transformed evolution of the brane z, and the parameters eu oc exp (fi — cj>). Notice that the latter 
satisfy <C 1 while the brane is evolving, and that the brane reverses direction whilst still strictly inside the orbifold interval. 



behaviour of v. Thirdly, the original SU(3) field v can already be made to reverse direction for particular choices of 
constants (see Eq. Hence, once mixed, the transformed z solution also inherits this bouncing behaviour. 

As before, this behaviour is subject to a' corrections. However, the strength of these corrections can always be 
adjusted such that, when the brane is bouncing, the corrections are extremely small and so under control. Of course, 
the corrections cannot be made arbitrarily small for all time, but they can always be made arbitrarily small over 
significant periods of time when the brane is moving. To achieve this one simply tunes Oi -C 1, which has the effect 
of setting (3 ^> in the vicinity of the bounce. 

These bouncing solutions richly extend the results of the previous sections. We now see that the effective supergravity 
action Eq. Q admits exact solutions where the brane evolves in a regime with £j<l, has small a 1 corrections, is 
strictly between the boundaries, and can also reverse direction mid-course. These effects were not at all obvious from 
the exactly integrable SU (2) and SU(3) Toda models, and yet can be generated by judiciously applying symmetries of 
the equations of motion. We also reiterate that no explicit potential was required to induce these effects; the reversal 
is a natural outcome of the non-linearly coupled cosmology. 



IX. PERTURBATIONS 



In the previous section we presented several new classes of cosmological solution to heterotic M-theory, and found that 
the four-dimensional scale- factor a = exp (a) always satisfies a ~ \t — iol 1 ^ 3 - Switching to conformal time r\ defined 
by di] = a~ 1 dt, this translates into a ~ 1 77 — 770 1 1 ^ 2 . This means that on the (+) branch we always have an expanding, 
decelerating universe, whilst on the (— ) branch we always have a contracting, inflationary universe. This behaviour 
is to be expected, for the cosmology we arc studying has no explicit potentials and so a remains unaffected by the 
other fields. We will now consider in more detail the inflationary epoch and the generation of perturbations on the 
(— ) branch. 

As with the familiar PBB scenarios, the inflationary period on the (— ) branch is characterised by a comoving Hubble 
length \d(hia)/d'r/\~ 1 = 2 1 77 — T70 1 that decreases as we take 77 — ► tjq and approach the Big Bang singularity from below. 
Consequently, a given comoving scale starting inside the Hubble radius as r) — > — 00 automatically becomes larger 
than the Hubble radius as r/ — > r?o- Therefore, on the (— ) branch one can produce super-horizon scale perturbations 
merely from kinetic-driven inflation, without the use of any potentials. This is considered an interesting alternative to 
conventional inflation on the (+) branch, since PBB scenarios do not require special choices of potential or slow-roll 
conditions. Given this, it is interesting to consider the perturbation spectra of our fields on the (— ) branch, and see 
whether there any useful scale-invariant modes. Not surprisingly, we will be able to utilise the techniques developed 
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in PBB cosmology to aid our calculations. We will also see that the factor of 3 in the kinetic terms for (3 and x, which 
complicated the classification of the scalar-field manifold (see Ref. ^(j) nas interesting consequences for the spectral 
indices of the fields. 

We will begin by considering perturbations around a special background where the axions x, a, v have been set to 
constants, and where the conserved quantity in Eq. I|15f) has been set to zero. In this case the brane z remains static 
for all time, and entirely decouples from the equations of motion of (3 and <j). The fields /3, 4> then exhibit standard 
rolling-radii behaviour with unrestricted parameters, and no transitions on the ellipse occur. Although this special 
"vacuum" situation does not incorporate any interesting brane displacements in the background, it proves to be a 
much simpler situation that can be solved analytically. Later, we will comment on perturbations around more general 
backgrounds, including the various Toda models and their symmetry transforms. In the meantime, we note that in 
the simple vacuum case the (3 and </> perturbations remain coupled to the metric perturbations, and produce adiabatic 
perturbations with the same steep n = 4 blue spectra that occurs in PBB cosmology [Tfilll7l|. In contrast, the fields 
z, a, Xi v with constant background values are decoupled from the metric perturbations, and produce isocurvature 
perturbations 6z, da, 5x, 5v with different spectra. 

The first-order, gauge-invariant perturbation equations for Sz, da, Sx, in conformal time are given by 

Sz" + (2a +[3' - 4>')5z' + k 2 Sz = 

5a" + 2(a! - <j)')5a' + k 2 Sa = -4qz((3' - <f>')5v' + Sqz X (3'Sz' (63) 

8x" + 2(a' - (3')5x' + k 2 S X = 
5v" + (2a' - (3' - 4>')8v' + k 2 5v = -2x/3'Sz'. (64) 

Here a / denotes a derivative with respect to rj, and k is the comoving wavenumbcr of the perturbation. In order 
to solve for these four isocurvature perturbations we will use techniques familiar from the PBB literature, with z 
replacing an axion. This involves making an appropriate conformal transformation on the metric into each axion's 
frame so as to eliminate the coupling to (3, 4>, and then solving the resulting perturbation equations in the usual manner 
(see Refs. EHHHl). 

However, before we do this we need to deal with the awkward source terms on the right-hand 
sides of Eq. (|63|l and Eq. 1)64(1 . The presence of the bulk- brane field z on the right-hand side of Eq. 1(63(1 , rather than 
a true axion, slightly complicates the situation as we cannot simply set z = as we can with x m Eq. 1(64(1 . Recall 
that our theory is only valid when z £ (0, 1). Instead, we must deal with the source terms by choosing appropriate 
combinations of perturbations: 5 A — Sa + 4qz5v and SB = 5i> — xSz . 

Following the calculations of Ref. we can now define a new metric for each field's frame by making a conformal 
transformation on the Einstein metric: (jj ^ v = fl'jg^. Our conformal factors flj are explicitly given by 

n 2 z =e^, n 2 A = e- 24 ', tf x =e- 2 P, 0| = e"^ (65) 

These conformal transformations lead to a different scale factor a,j = VLjd in each frame, depending on each fields 
coupling to (3, 4>. As we are considering static axions and bulk brane, (3, 4> behave as simple rolling-radii fields with 
fixed parameters that lie at one point on the ellipse, Eq. 1(22(1 . for all time. 2 Explicitly, in conformal time they satisfy 

i 3 3 

a = a* M 5 , P= 2^/3 ln M +A) , <t> = -jP<l> In \r)\ + (f>o (66) 

where a* is a constant, and we have conveniently set 770 = 0. Using these background solutions one can show that 

1 3 3 3 3 

a,j = a*j|?7|2+ r i with r z = - (pp - p<j,) , r A = --p 4 ,, r x = --pp, r B = -^(Pp +P<j>) (67) 

where the a*j are a set of constants. In these new frames we then find that the perturbation equations can be recast 
in the form 

5x'< + 2a' j 6x' j + k 2 5xj = 



1 Unlike z, we are always free to set x = 0. As such, we can choose SB = 5v. 

2 Consequently, we will now drop the subscripts i, f that label the initial and final rolling-radii powers, as the fields /3,<f> remain in the 
same rolling-radii states for all time. Also, note that we can pick {pj3,p^,) from anywhere on the ellipse, irrespective of the time branch. 
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where 5xj = (5z,5A,5x,5B) and a'j = a'j/dj is the Hubble rate in each conformal frame. The solution for our 
isocurvature perturbations, after normalising at early time, is then given by (see Ref. [T3 p 



5x4 



nijk 



cxp 



(1 + %I) 



a, l r il V " 



Here the rrij are given by m z = 2q, ttia = 4, m x = 12, mj = 8g, and Hj 1 ^ is the Hankcl function of the first kind and 
order J. Defining the power spectrum P$ x and its spectral index n$ x for a general perturbation Sx as 

P fe = — |fo| and n^-^—- 

we find the spectral index for each of the isocurvature perturbations is given by 

nsxi = 4 - 2|r,| 



Looking at the definitions of the \rj \ given in Eq. (|67|l . we sec how the spectral indices are dependent on the coupling 
of z and the axions to <f>,j3 and consequently their expansion powers, p^ and pp. Inserting the specific couplings for 
each field and considering the range of background solutions yields 



n 5A = 4-3|iv|: e[4-2v / 3,4] ~ [0.54,4], 
3 

nsz =4- -\(pp -P4,)\, ri4 x = 4 - 3|p/3| , n SB 



-|b/3 +P<t>)\ 



e [2,4] 



Thus wc find that our perturbation 5 A has the classic axion perturbation spectrum familiar from PBB calculations, 
and can provide a scale-invariant spectrum. In contrast, the bulk brane and other axion perturbations cannot provide 
a scale-invariant spectrum, a result similar to the one obtained in Ref. [l3T | . 

One can also write the spectral indices as a function of a single variable by using the ellipse constraint, Eq. I122II . This 
reveals 



n Sz = 4 



±J4-3pd-3p, 



, 71^ = 4-31 



n Sx = 4 - ^4-3^ , n§B = 4 - 



±,/4-3*H +3p* 



where p G [-2/V3, 2/V3] for the vacuum case. One should remember that the choice of ± sign must be consistently 
applied across all the spectral indices, and that both signs are always valid choices (as we do not have to satisfy 
Pfj ~ P(f> > in the vacuum case). 

If one is familiar with PBB calculations the above result may be surprising, as axion perturbations derived from 
actions very similar to ours will usually all have spectral indices in the range [4 — 2-\/3, 4]. (See, for example, the 
variety of dilaton-moduli-axion systems discussed in Refs. 0,0,0,^3). This change is a direct consequence of the 
coupling of our fields to f3, which unlike in PBB cosmology has a factor of 3 in its kinetic term. This then affects the 
range of pp through the ellipse condition. One cannot change this result by rescaling /3's kinetic term as this rescales 
/3's coupling to the fields and moves the effect into the rj definitions. This then leaves 5 A as the single perturbation 
capable of producing a scale-invariant spectrum. 

So far we have only been considering the "vacuum" solutions where z and the axions remain constant. However, 
generalising these solutions to the case with moving brane and axions remains an open question, due to the sheer 
complexity of the solutions considered. One can begin by truncating off the axions and considering the perturbations 
of the 4>,/3,z SU(2) action of Section llVl In this case one can use an SL(2,M.) symmetry of the truncated action 3 
to solve for the perturbation Sz around a mewmg-brane SU(2) background, by applying the <SX(2,R) symmetry to 
perturbations 5(j>, 5(3, 5z around a static-bieme SU(2) background. One then finds that this "rotated" 5z isocurvature 
perturbation retains the spectrum n 6 [2,4] (see Ref. [l3|)- However, the effect of this rotation on the remaining 
axionic equations leaves a non-trivial calculation. 



3 This symmetry is not related to the SL(2, R) generalised T-duality we have discussed in this paper, and does not remain when considering 
the full, untruncated action. 
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As a result, we can only conjecture that a scale-invariant mode persists in perturbations around the Toda model 
backgrounds and their symmetry transforms. However, it is certainly true that all of the solutions we have considered 
will asymptotically approach the "vacuum" scenario. Moreover, when one applies the constraints on p^ in the various 
classes of solutions one finds that, in a least one of the asymptotic limits, the rolling-radii regime which leads to SA 
producing a scale-invariant spectrum is accessible. Hence, we can always generate a scale-invariant mode in one of 
the asymptotic limits, even if we cannot determine whether such a mode can also be generated at intermediate times. 



X. CONCLUSION 



We have presented several new classes of cosmological solution to the four-dimensional effective supergravity 
description of heterotic M-theory. This theory contain seven fields: the four-dimensional scale- factor a, the modulus (3 
measuring the separation of the orbifold planes, the axion x related to the graviphoton field, the dilaton 4> measuring 
the average Calabi-Yau volume, the axion a related to the bulk three-form, the field z locating the position of the M5 
branc, and the axion v representing the self-dual two-form on the branc worldvolumc. To linear-order in the moduli- 
dependent parameters (fc = 1, 2), all fields except a can be described by the following scalar-field Lagrangian 



C = 



\d(3) 2 + 3e- 2f} (d X ) 2 + ](d(j>) 2 + -e^ (8a + Aqzdv) 2 + \qe^^{dz) 2 + 2qe- fl -*{dv - X dz) 2 



We have attempted to identify as many exact solutions to this system as possible, by identifying special constraints on 
the fields that simplify the analysis. The only previously known solution to this Lagrangian, as described in Ref. Q, 
is found when C is consistently truncated to the form 

Csu® = \{d(3) 2 + \W) 2 + \qeP-+{dz? 

The fields /3, <fi, z then form an exactly-solvable SU (2) Toda model, with the brane z undergoing single displacements. 
In this paper we have identified three new solutions in addition to this ST/ (2) Toda solution. The first new solution 
was found by consistently truncating C to the different form 

C SU{3) = \{dP) 2 + 3e- 2 W + \{d<P) 2 + \qe^{dz) 2 

using the conditions dv — xdz = da + Aqzdv = 0. By switching off these two terms, one finds that the reduced set 
of fields /3, x, 4>i z span an integrable St/(3) Toda model, and can be solved for exactly. This SU(3) model allows for 
double displacements of the brane z, and the SC/(3) solutions can always be made reliable with e^Cl over a certain 
period of time during this double displacement. However, the brane must leave the compact space in any solution 
that has reliable £k <C 1 regime at some point. 

Next, we applied to the SU(2) and SC/(3) models the symmetry transformations derived and discussed in our 
companion Ref. [Tpj . This enabled us to derive two new and distinct cosmological solutions. The properties of 
these new solutions were then discussed at length, and it was found that the reliability of the solutions had been 
radically affected. This is ultimately due to the SL(2,M) subgroup of symmetries, which acts as a "generalised" set 
of T-duality transformations given by 

T , = aT_ib z , = _^_ S < = S _2^L_ with aj& ,c,d€R , ad - be = 1 (68) 
icT + d icT + d icT + 1 

This is not ordinary T-duality, as all three complex fields T, S, Z arc affected. Using these symmetries, it was then 
found in the SU (2)-transformcd solutions that the brane can undergo a single displacement entirely within the orbifold 
interval with -C 1 throughout. It was also found in the S£/(3)-transformed solutions that the brane field z can 
undergo two successive displacements of opposite sign and so reverse direction. The specific conditions under which 
this reversal occurs are as follows. Firstly, set da + Aqzdv = so that the axion a is decoupled from the other fields. 
Then the scalar-field lagrangian L reduces to the simpler form 

£' = lid/3) 2 + 3e- 2l3 (d X ) 2 + \{dcfy 2 + ^qe^(8zf + 2qe~^{dv - X dz) 2 

One then proceeds by setting dv — xdz = and solving the system as an SC/(3) Toda model, but then restoring 
the dv — xdz term to a general, non-zero value by applying an SL(2, R) symmetry. In particular, the fields z,v 
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transform as a doublet under SL(2,M), and so we can solve the system with general dv — \dz by "rotating" from 
a solution where it is zero. As a consequence, the equations of motion arising from the reduced lagrangian £ have 
been completely solved in this paper. Further, by tuning the sign and magnitude of the -~x^ z contribution generated 
by the symmetry application, one can modify the overall velocity of the brane so that it comes to rest and reverses 
direction. This is a particularly interesting feature arising from the coupling with Xi whose presence in the kinetic 
term dv — \dz is due to the need for a gauge-covariant derivative in five dimensions. This reversing behaviour, which 
we have occasionally called a "bouncing" solution, can also be made to occur entirely within the orbifold interval with 
£fe <C 1 throughout. 

As such, all of the transformed solutions demonstrate a rich new variety of M5 brane behaviours, and new, trustworthy 
regions of solution space emerge that had not previously been identified. In particular, we conjecture that reversing 
brane solutions will exist in other corners of string theory beyond heterotic M-theory. One can pin down reasonably 
clear "minimum conditions" for this reversal to occur, as follows. Firstly, at least one modulus should be active, such 
as the dilaton </>, whose coupling to the brane kinetic term will induce the brane z to undergo a single displacement. 
Secondly there should also be an active combination proportional to a cross-coupling between an axion field and dz. 
This second combination can then be adjusted so that the brane turns around at some point during its motion. 

As an interesting corollary, we then considered the isocurvature perturbation spectra produced by the model in an 
inflationary contracting (PBB) phase. In the "vacuum" case we found that one of the isocurvature modes - the one 
associated with the axions a and v — is able to produce a scale-invariant spectrum. Furthermore, we found that all 
of the solutions considered will asymptotically approach this vacuum case in at least one asymptotic limit, and so 
a scale-invariant perturbation spectrum can always be generated asymptotically when perturbing around any of the 
solutions we have studied. However, the detailed structure of the perturbation spectrum at intermediate times has not 
yet been computed in its full generality, and it would be interesting to study this problem in greater depth, perhaps 
in a manner analogous to the numerical approach developed in Refs. |l9ll20l|. 

Finally, we note that the methodologies employed in this paper have much wider applicability. For example, the Toda 
model solution method, as extensively detailed in Refs. jg, Q , is not restricted to scalar-field systems arising from 
heterotic M-theory, and could be readily utilised in other areas of string theory. Likewise, it is equally plausible that 
other braneworld Kahler metrics may possess useful symmetry groups, which can be used to transform subsystems 
of fields into new patterns of behaviour. In light of this, it would be interesting to clarify the origin of the special 
SX(2,R) symmetry group that we have found, and understand the general conditions under which reversing brane 
behaviour occurs in string and M-theory. 
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Appendix A 



In this Appendix we present certain additional details of the SU(3) Toda model derivation. We do this because 
the derivation is rather complicated, particularly the manner in which one must change time-gauges and judiciously 
redefine constants. 

To begin with, we know that the vectors qi,q2 are proportional to the two simple root- vectors of SU(3). Utilising 
this fact, we can choose a basis for the space (a,f3,<fi) that is adapted to the underlying SU(3) symmetry, and so 
consists of vectors e ,ei,e 2 satisfying 

3 3 

<e ,e >=-l , < e ,ei > = < e ,e 2 > = , ei = -qi , e 2 = -q 2 (69) 
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A choice of basis compatible with these conditions is given by 

3,„ 3 



eo 



(73, 0,0) , ex = -(0,-1,1) , e 2 = -(0,2,0) (70) 



We now write the covariant vector Gat as the following sum 

2 

Ga = ^ Pl (T)e, (71) 

i=0 

and insert this time-dependent expansion into the equations of motion to find the evolution of the "modes" pi. 
Choosing the convenient gauge n — 3a (or E = 1) one finds 

Po = (72) 
p\ + ^ufe 2 "-* = (73) 

92 + \u 2 2 e 2p2 ~ pl = (74) 

-Po + \{pI- PiP2 +pt)+2U = (75) 
The general solution to these equations is now easy to come by, and takes the form 

Po = -ko{T-T ) (76) 
Pi = -ln fl i(r) (77) 
p 2 = -ln 52 (r) (78) 

where fco,7o are constants. The functions gi,.g 2 are given by a sum over the collection of weight vectors 
Ai = {(0,-1), (-1,1), (1,0)}, A 2 = {(-1,0), (1,-1), (0,1)} of the fundamental 3 and 3 representations of SU(3). 
Concretely, if we define the matrix of vectors 

(ArWtO,-!) (-1,1) (1,0) \ 
I A 2 ) I (-1,0) (1,-1) (0,1) ) 1 j 
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then for i = 1, 2 the functions are given by 



3 

■9* = X] a y exp ' ( kr - r )] ( 80 ) 

J'=l 



where the positive constants are 

4m 2 / 2fci - fc 2 



ail = — I" — 

,2 



4uj fki + k 2 



fl 12 = — ^ p 

,2 



4wf / 2fc 2 - fci 



a 13 = — I p 





(Ih-i-kx' 




\ P 




(ki + fc 2 ^ 








^2fci - fc 2 







and P = (2fci — fc 2 )(2fc 2 — fci)(fci + fc 2 ). The constant vector r = (ri, r 2 ) is a set of arbitrary time-shifts. The constant 
vector k = (fci, fc 2 ) is restricted to the open Weyl chamber, which means it is forced to have positive scalar product 
with the two simple root-vectors as follows 



(fci.fe) • (2,-1) = 2fa -hi > 
(fci,*3)-(-l,2) = 2fe-*i >0 



(81) 



These two conditions guarantee that gi,g 2 > so that the logarithms in pi,p 2 are always well-defined. Lastly, we 
must also impose the Friedmann constraint 







(82) 



Using the information above, one can arrive at an explicit solution for the fields a, j3, 4> and the two additional fields 
z, x that were integrated out. Recall that 



and that 



Then wc find 



a 




3i (r) = a n e- k2T+T2 + a 12 e^- k ^ T -^- T ^ + a 13 e klT - Tl 
92{t) = a 2 ie- fclT+ri + a 22 e^- k ^ T -^- T ^ + a 23 e k * T - T > 



(83) 



(84) 



Ot = -y= k) (T - T ) 



P= -In 

-In 

3 



one 



— k2T + T2 



+ a 12 e {k2 - kl)T - {T2 - Tl) + a 13 e klT - 



(12 



ie -kir+r 1 + a22e(i:i -fc 2 ) T -( Tl - r2 ) + ^^fc, 



zo = C : 



hi 

2 

9i 

"'9i 

<)2 



(85) 



a lie - k2T+T2 + a 12 e( k2 - kl)T -<> T2 - Tl) + a 13 e klT - Tl 

k 2 a lie - k2T+T2 + (k 2 - fci)ai 2 e( fc2 - fel ' T ~( T2 - Tl ) + k iai3 e klT - Tl 



C z 



X-Xq = C x — = C x 

92 



aiie- fe2T+r 2 + a 12 e( fc 2- fe i) T -( T 2- T i) + ai 3 e fe i r - r i 
-k 1 a 21 e- klT+Tl + (fci - k 2 )a 22 e^ kl - k2 ^ T -^ 1 - T2 ^k 2 a 2S ,e k2T - T2 
a 2 ie~ klT+T1 + a 22 e( fcl " fc2 ) T -( T1 - T2 ) + a 23 e k2T ~ T2 
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where zq, xo are constants of integration, and 



C z 



( 9 



Vou\ 



\8qu 2 

We can then deduce the forms of a, v that are compatible with the ancillary conditions Eq. i|26|) 

"fc 2 (2fc! - fc 2 )a 13 e felT - ri + k x {2k 2 - k 1 )a 23 e~ k ^+^- 



v - vq = xa(z - z Q ) + C Z C X 



cine 



+ T2 + ai 2 e( fe 2- fc i) T -( T 2— n) + a 13 e klT - Tl 



a - (Tq 



-2q 



2z v + xo(z - z ) 2 + / [(z - z Q ) 2 ]' (x - Xo)dr 



where vq,ctq arc two further constants of integration. Notice that the integral in the a solution is not elementary, 
and so cannot be written as a finite (potentially nested) sequence of logs, exponentials and rational functions of r. 
However, it can sometimes be analytically integrated for fixed choices of the arbitrary constants, and in any event 
has a sensible definite integral between fixed r limits. In particular, one can compute the a behaviour numerically for 
any given set of starting conditions. 

It is useful to understand the asymptotic limits, in order to transform these solutions to the proper time-gauge n = 0. 
One finds that 

e P „ e (2 kl - k2 )T/2 and e ^ e 3fe 2 r/2 ag T ^_ QO 
e /3 ^ e ( fcl -2fe 2 W2 and e 0^ e -3fc ir /2 ag r ^ +QO 



Since 2fci — k 2 > and k\ — 2k 2 < we see that the orbifold radius /? always goes from a state of expansion at early 
time to a state of contraction at late time. The same is true of the modulus cf) that measures the orbifold-averaged 
Calabi-Yau volume. These two fields will then have some complicated intermediate transition(s) that smoothly link 
these extreme limits. On the other hand, the fields z and x always asymptote to constants in the limits, although 
these constants are generally different. They too will undergo some intermediate "displacement" consistent with the 
different constant field values at early and late time. 

We now change from logarithmic time r to proper time t. Since the logarithmic-time gauge is given by n = 3a, we 
can find the relation to proper time by integrating the defining relation 

dt = e n{T Ur = e 3Q(T) dr = e^^^W (86) 

This gives 



In 



V3ko(t-t ) 





t-to 




3 





(87) 



where to is a finite integration constant. This leads to two disconnected time branches corresponding to the choices 
t — to > 0,Xb > and t — to < 0,Tb < 0, both of which lead to a well-defined positive argument for the logarithm. 
The regime t — to > 0,Xb > will be referred to as the "positive-time" or simply (+) branch, whilst the sector 
t — to < 0, To < will be dubbed the "negative-time" or simply (— ) branch. The physics in the time interval 
t G (-co, +oo) is mapped to these two regions in the following way. The early time r — * — oo regime with expanding 
/3, 4> corresponds to t — to —* on the (+) branch and t — to — > — oo on the (— ) branch, while the late time regime 
t — ► oo corresponds to t — to — > oo on the (+) branch and t — to — > on the (— ) branch. It should be noted these two 
time branches in t are physically separated by an unavoidable curvature singularity at t = to, despite the fact that in 
the r-gauge we had only one physical region. After this gauge change, typical terms in p\ , p 2 will then scale as 



t-to 


K 


t-to 




To 




To 





respectively, where A • = (k[ , k 2 



k[, —k' 2 ) and 



k[ =k x To , k' 2 = k 2 T Q 



25 



Note that k[,k' 2 can be of either sign depending upon the sign of To and so the choice of branch. On the positive 
branch we have k[, k' 2 > so that k[ = \k[\, k' 2 = \k' 2 \. In this gauge the fields {3, 4> scale asymptotically at early-time 
as 



where we have defined the two constants 



1 



(2|*i|-|^|) 



\k'\ 



The Friedmann constraint Eq. I|82|) then reduces to the familiar ellipse condition 



Moreover, the Weyl chamber constraints 2\k[\ — \k' 2 \ > 0, 2\k' 2 \ — \k[\ > translate into pp^ > 0, 5 < 0. So the 
positive branch is associated to 6 < with the further additional constraint that 8p = —2ppj < 0. Conversely, on the 
negative time branch we find the opposite results, since we follow the ellipse trajectories backward. Hence, we find 
that 5, Sp > 0. 

Since all the fields involved are scalars, we are now free to substitute r in terms of t in all the solutions. These will 
then be the n = forms for the solutions. To make these solutions look "nice" , however, one must carefully redefine 
certain constants. If one defines new timcscalcs T, Tp via 



— exp [-(ki - 2k 2 )r Q + (n - 2r 2 )] 
an 



— exp \(2ki - k 2 )r - (2t\ - r 2 ) 
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then the resulting solutions look relatively simple. They can be made even simpler by defining the ubiquitous fractional 
combinations 



013021 



023011 



012022 



012&22 



Finally, we have decided to write the solution for {3 so that it scales with respect to T at early time rather than the 
natural choice Tp, and so the ratio T/TJg has been absorbed into the additive offset (3q. This brings the (3 solution 
closer to the old ST/ (2) form and simplifies the vector notation, but at the expense of introducing the ratio T/Tp into 
the constraints. 



